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Abstract 

^Prom the analyticity properties of the equation governing infinitesimal per- 
turbations, it is shown that all stability properties of spatially extended ID 
systems can be derived from a single function that we call entropy potential 
since it gives directly the Kolmogorov-Sinai entropy density. Such a function 
allows determining also Lyapunov spectra in reference frames where time- 
like and space-like axes point in general directions in the space-time plane. 
The existence of an entropy potential implies that the integrated density of 
positive exponents is independent of the reference frame. 
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I. INTRODUCTION 



The first part of tliis work (liereafter referred to as LPT) was devoted to the definition 
and discussion of the properties of temporal (TLS) and spatial (SLS) Lyapunov spectra of 
ID extended dynamical systems. In this second part we first show how the two approaches 
are mutually related by proving, in some simple cases, and conjecturing, in general, that 
all stability properties can be derived from a single observable: the entropy potential, which 
is a function of two independent variables, the spatial and the temporal growth rates /i. A, 
respectively. Legendre transforms represent the right tool to achieve a complete description 
of linear stability properties in the space-time plane. In fact, we find that equivalent descrip- 
tions can be obtained by choosing any pair of independent variables in the set {n^, tza, /i. A}, 
where and ux are the integrated densities of spatial, resp. temporal, Lyapunov exponents. 
The corresponding potentials are connected via suitable Legendre transformations involving 
pairs of conjugated variables. 

A further type of connection between spatial and temporal Lyapunov exponents discussed 
in this paper is found in connection with the evolution of perturbations along generic "world- 
lines" in the space-time, i.e. along directions other than the space and time axes considered 
in LPT. The extension of the usual definition of Lyapunov exponents to this more general 
class of frames, already discussed in , is rather appropriate for characterizing patterns with 
some anisotropy. Here, we show that this seemingly more general class of spatiotemporal 
exponents can be derived from the knowledge of spatial and temporal Lyapunov spectra, 
which thus confirm to contain all the relevant information. 

For the sake of completeness, finally, we recall the last class of exponents introduced to 
describe convectively unstable states, comoving Lyapunov exponents [Q, and their relation- 
ship again with SLS and TLS. In particular, we discuss the structure of the spectra in a 
simple case of a stationary random state. 

Let us now briefly introduce the notations with reference to some specific models. Spa- 
tiotemporal chaos and instabilities in extended systems have been widely studied with the 
aid of models of reaction-diffusion processes, whose general ID form is of the type [^,|| 



with the state variable y{x,t) defined on the domain [0, L] (periodic boundary conditions 
y(0,t) = y{L,t) are generally assumed). The nonlinear function F accounts for the lo- 
cal reaction dynamics, while the diffusion matrix D represents the strength of the spatial 
coupling. The introduction of coupled map lattices (CML) has been of great help for under- 
standing the statistical properties of spatio-temporal chaos, especially by means of numerical 
simulations. 

In its standard form a CML dynamics reads as 



where i, n being the space, resp. time, indices labelling each variable of a lattice of length 
L (with periodic boundary conditions ?/^+^ = y^), and e gauges the diffusion strength. The 
function /, mapping a given interval / of the real axis onto itself, simulates a local nonlinear 
reaction process. 



dty = F(y) + mly 



(1) 




(2) 
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A generalization of model (0) has been proposed |^ to mimic ID open-flow systems, 
namely 

= /((!- e)yl, + e[{l- a)yl' + ayi,+']) . (3) 

The parameter a (bounded between and 1) accounts for the possibility of an asjTiimetric 
coupling, corresponding to first order derivatives in the continuum limit. 

The present paper is organized as follows. In Section II we introduce the entropy poten- 
tial and derive its explicit expression in some simple cases. Sec. Ill is devoted to Lyapunov 
analysis in tilted reference frames, while Sec. IV deals with the relationships between tem- 
poral, spatial and spatiotemporal exponents. Comoving exponents are reviewed in Sec. V 
within the framework introduced in this paper. Some conclusive remarks are finally reported 
in Sec. VI. 



II. ENTROPY POTENTIAL 

The simplest context, where a discussion on the entropy potential can be set in, is 
provided by the linear diffusion equation for the field u{x,t) 

dtu = 7M + Ddlu , (4) 

which can be interpreted as the linearization of (the scalar version of) Eq. (|l|) around a uni- 
form stationary solution y{x,t) = const. The linear stability analysis amounts to assuming 
a perturbation of the form 



U[X 



t) ~ exp (^fix + Xtj , (5) 



where X = X + iu and jl = ^ + ik are complex numbers the real parts of which denote 
temporal and spatial Lyapunov exponents, respectively, while the imaginary parts represent 
the integrated densities of spatial {u) and temporal {k) exponents. Substituting Eq. (H) in 
Eq. (^ we obtain 

A = 7 + D/i^ (6) 
By separating real and imaginary parts, we get 

X{uj,k) = -i+iuo^ -AD^k^) /{ADk^) (7) 



fx{iu,k) = uj/{2Dk) . 

As already discussed in LPT, uj and k play the same role as the integrated densities ra^ and 
nx and can be explicitely obtained by inverting Eqs. (H), 



= = - (8) 
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The minus sign in the definition of is just a matter of convention: we adopt this choice 
for consistency reasons with LPT. 

The above sets of equations (|^||) stress, in a particular instance, the general observation 
reported in LPT that either the pair {n^^nx) or (/i. A) suffices to identify a given perturba- 
tion, the remaining two variables being determined from the Lyapunov spectra. However, 
any two items in the set {yU, \,n^j_,n\} can be chosen to be the independent variables. The 
above two choices are preferable for symmetry reasons; however, the pairs (/x, n\) and (A, n^) 
turn out to be the best ones for the identification of a single function, the entropy potential, 
which determines all stability properties. 

In fact, as it is clear from Eq. (^, we can condense the two real functions needed for a 
complete characerization of the stability properties into a single analytic complex expression. 
Now, the mere circumstance that A(/x, n\) and n^(/i, nx) are the real and imaginary parts of 
the analytic function A*(/i) has an immediate and important consequence:Q Cauchy-Riemann 
conditions are satisfied and it is possible to write A and as partial derivatives of the same 
real function, 



a* 

dnx 
a* 
9/i 



A (9) 



where \l/ is the imaginary part of the formal integral \l/ of A with respect to jl. Equivalently, 
one might call into play the real part of \E', as it is known that the latter contains the same 
amount of information. 

In the case under investigation, we find 

^(/i, nx) =nx{i + D^^) - ^nl , (10) 

which, together with Eq. (|^), provides a complete characterization of the system. 

Another, less trivial, example where the linearized problem leads to an analytic function 
for the eigenvalues is the the ID complex Ginzburg-Landau equation 



dtA = (1 + ic{)dlA + A - (1 - zc3)A|A|2 , (11) 

where A{x, t) is a complex field and Ci and C3 are real positive parameters. The stability 
of the "phase winding" solutions A{x,t) = AQexp^i^ux — uqI)), with Aq = \J\ — v'^ and 
tuo = — C3 + (ci + Cs)//^, are ruled by the following equation for the (complex) perturbation 

m(x, t) 

dtU = {\^%Cx){dlu^2id^u)-{\-ic^){\-v^)[u^u*) , (12) 

together with its complex conjugate for u* considered as an independent variable. The 
eigenvalue problem is solved assuming again 



^The reference to the complex conjugate variable again follows from the convention adopted for 
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u{x,t) = Uoexp ynx + Xtj ; u* {x , t) = exp yfxx + Xtj , (13) 

and equating to zero the determinant of the resulting hnear system, to get the analytic 
(implicit) relation between A and ft 

{X + l-u^-jl^ + 2ciujiy + (ci/i' + 2z//i + (1 - z/2)c3)' = (1 + cl){l - z/2)2, (14) 

which is analogous to Eq. (||) 

On the basis of the examples discussed here and in the Appendix, one can convince 
himself that the analyticity property of the eigenvalue equation appears to be very general. 
Periodicity in time simply leads to multiply several r.h.s.'s all depending on ft, while period- 
icity in space requires distinguishing between different sites on the lattice. In the latter case, 
the equivalent of Eq. ( |A7| ) is obtained by equating to zero a suitable determinant, where 
the only variable is jl. There is no reason to expect that different conditions should hold in 
aperiodic regimes. 

The above approach is, in some sense, a generalization of dispersion relations which are 
normally introduced for the characterization of elliptic equations. In that case, the only 
acceptable linear solutions are propagating plane waves, that is A = /i = for (almost) 
all wavenumbers. From our point of view, this implies a strong simplification since the 
mutual relationships among {nx,n^, X, ^} reduce to the link between spatial and temporal 
wavenumbers. Moreover, it is obvious that, because of the degeneracy, not all representations 
are equivalent (in particular, the (/i. A) plane is totally useless). 

We must stress that the methodology that we are trying to develop in these two papers 
applies to general systems where propagation coexists with amplification (or damping). This 
is by no means a limitation, as all models introduced for the characterization of space-time 
chaos are in this class. 

However, the most serious obstacle to a rigorous proof of the general validity of Eq. (^ 
is represented by the identification of the integrated densities n^, nx with the wavenumbers 
CO and k, respectively. In the presence of spatial disorder, the Lyapunov vectors are no 
longer Fourier modes: one can at most determine an average wavenumber by counting the 
number of nodes in the eigenfunctions. This is not a problem in the absence of temporal 
disorder, when the node theorem applies 0. However, in more general cases, it is no longer 
possible to speak of eigenfunctions and we are not aware of any generalization to overcome 
the difficulty. For such a reason, we have performed some direct numerical check to verify 
the correctness of our conjectures. 

Before discussing numerical simulations, let us come back to the problem of the represen- 
tation. In the above part, we have seen that the choice of the pair of independent variables 
(/i, n\) was very fruitful for the identification of a potential. However, the asymmetry of such 
a choice calls for transferring the above result in either representation proposed in LPT. This 
step can be easily done with the help of Legendre transforms. We discuss the transformation 
to the plane (/i. A), any other transformation being a straightforward generalization of the 
same procedure. 

/^From the first of Eq. (^, we see that A and nx can, indeed, be considered as conjugate 
variables in a Legendre transform involving The conjugate potential is naturally 

$ = AriA - ^ . (15) 
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It is easily seen that in the new representation, the following relations hold 

= riA (16) 

Accordingly, the potential $ is the appropriate function which allows determining the two 
integrated densities in the symmetric representation (A,/i). We call $ the entropy potential 
since it coincides with the Kolmogorov-Sinai entropy density along a suitable line (see Sec. 
III). 

In Fig. |I], we present a numerical reconstruction of $(/i. A) in terms of its contour levels for 
the homogenoeus CML, namely model (0) with f{x) = rx(mod 1). The entropy potential 
is obviously known up to an additive arbitrary constant that we have fixed by imposing 
that the value attained on the upper border is equal to zero. The potential increaeses 
monotonously from top to bottom. It is clear that outside the allowed region delimited by 
the solid curves, $ is a linear function of n and A. 

The structure of the potential does not substantially change for more general CMLs. 
We have tested Eq. ([T6|) for a lattice of logistic maps (/(x) = 4x(l — x)) with e = 1/3 by 
integrating along two different paths in the (/x. A) plane (see Table 1). The difference is so 
small that we can confirm that the relations are valid, within the numerical error. 



III. SPATIOTEMPORAL EXPONENTS 

In the perspective of a complete characterization of space-time chaos, one should consider 
the possibility of viewing a generic pattern as being generated along directions other than 
time and space axes. In fact, once a pattern is given, any direction can, a priori, be considered 
as an appropriate "time" axis. Accordingly, questions can be addressed about the statistical 
properties of the pattern when viewed in that way. 



A. Definitions 

For the moment, we assume that the pattern is continuous both along space and time 
directions; we shall discuss later how the definitions can be extended to CML models. Let 
us consider a given spatiotemporal configuration of the field y{x,t), generated, say, by inte- 
grating Eq. (p. When arbitrary directions are considered in the (x, t) plane, the coordinates 
must be properly scaled in order to force them to have the same dimension. We choose to 
multiply the time variable by c, where c is a suitable constant with the dimension of a ve- 
locity. Moreover, let "i? denote the rotation angle of the tilted frame {x', ct') with respect to 
the initial one {x,ct), adopting the convention that positive angles correspond to clockwise 
rotations. Sometimes, it will be more convenient to identify the new frame by referring 
to the velocity v = ctan^}. The limit cases v = {-^ = 0) and v = +oo ("i? = vr/2) cor- 
respond to purely temporal and purely spatial propagations, respectively. The coordinate 
transformation reads as 

ct' = (3(^ct + ^x^ (17) 
x' = (3 {-vt + x) , 
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where (3 = 1/y (1 + f^/c^). The physical meaning of v is transparent: it can be interpreted 
as the velocity in the old frame of a point stationary in the tilted frame (constant x'). 

The new field y{x',ct') can be thought of as being the result of the integration of the 
model derived from the original one after the change of variables Although it is not 

obvious whether the invariant measure in the initial frame is still attracting in the new frame 
(see Ref. [0] for a discussion of this point), one can anyhow study the stability properties by 
linearizing and defining the Lyapunov exponents in the usual way. 

In CML models, the discreteness of both the space and the time lattice leads to some 
difficulties in the practical construction of tilted frames. In fact, only rational values of the 
velocity v can be realized in finite lattices (in this case, it is natural to assume that the 
lattice spacing is the "same" along the spatial and the temporal directions and, accordingly, 
to set c = 1). Moreover, writing the explicit expression of the model requires introducing 
different site types. For this reason, we discuss in the following the simplest nontrivial case 
f = 1/2, the generalization to other rational velocities being conceptually straightforward. 

A generic spatial configuration in the tilted frame is defined by sites of the spatiotemporal 
lattice (^,n) connected by alternating horizontal (as in the usual case) and diagonal bonds 
(see Fig. 0). By suitably adjusting the relative fraction of the two types of links all rotations 
between and 7r/4 can be reproduced. The explicit expression of the updating rule requires 
a proper numbering of the consecutive sites. Moreover, as seen in Fig. 0, it involves the 
"memory" of two previous states. 

Finally, an exact implementation of the mapping rule requires acausal boundary condi- 
tions, since the knowledge of future (in the original frame) states is required (this is a 
general problem occurring also in the continuous case). As we are interested in the thermo- 
dynamic limit, we bypass the problem by choosing periodic boundary conditions. Such a 
choice has been shown not to affect the bulk properties of the dynamical evolution [0. 

In the updating procedure, two different cases are recognized: the variable y is either 
determined from the past values in the neighbouring sites, or it requires the newly updated 
I/- value on the right neighbour (see Fig. For v = 1/2, this can be done by simply 
distinguishing between even and odd sites. 



^nVi = /((!- ^)Xl' + I [y^'-' + X':^' ] (18) 



where i = 1, . . . , L/2 (L is assumed to be even for simplicity), while 

YUi^Xi , (19) 

are additional variables introduced to account for the dependence at time n — 1. Taking into 
account that X'^^^ can be determined from the X and Y variables at time n, the mapping 
can be finally expressed in the usual synchronous form (X^, Y^) (-^^^+1, ^n+i)' but with an 
asymmetric spatial coupling with next and next-to-next nearest neighbours. The Lyapunov 
exponents r/j can now be computed with the usual technique 0. 

In analogy with the original model, we expect again that, in the limit of infinitely ex- 
tended systems, the set of exponents rij{v) will converge to an asymptotic form. 
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(20) 



where is the corresponding integrated density. We will refer to this function as the 
spatiotemporal Lyapunov spectrum (STLS). In the limit cases f = 0, +00 {'d = 0, vr/2), the 
STLS reduces to the standard temporal and spatial spectrum, respectively. 

The recursive scheme (|^) implies an increase of the phase-space dimension by a factor 
(1+1/2) (in general 1 + v). Actually, as we will argue, these new degrees of freedom are 
not physically relevant. However, for consistency reasons with the original rescaling of the 
spatial variable, we choose to normalize the spatiotemporal density between and 1 + f (the 
time units are, instead, left unchanged by the above construction). 

The generalization to asymmetric maps (^ is straightforward: it removes the degeneracy 
V —>■ —V. Numerical results for logistic maps, indicate that the dependence of the positive 
exponents on the velocity is quite weak in the fully symmetric case a = 1/2 (for instance, the 
maximum exponent exhibits a 20% variation in the whole v range), while it is remarkable 
for asymmetric couplings. In every case, the negative part of the spectrum sharply changes 
with the velocity. This is consistent with the results obtained for delayed maps in Ref. [Q]. 



B. Representation in the {fi, A) plane 

Spatiotemporal exponents can be put in relation with fi and A by rewriting the general 
expression for a perturbation in a frame rotated by an angle ^, 

exp(/ix + At) = exp(/iV + X't') . (21) 

Such an equation induces a rotation of the same angle in the (c/i. A) variables, 

A' = /5 (A + vjj) (22) 
fi' = P {-{v/c^)X + fi) . 

The above equations allow studying the stability with respect to generic perturbations with 
an exponential profile along x'. For simplicity, we shall consider only uniform perturbations, 

exp{fj,x + At) = exp^rjt') , (23) 

where the growth rate rj denotes the spatiotemporal exponent. Notice that we have changed 
notations from A' to rj, to understand that the condition /i' = is fulfilled. From the second 
of Eq. (^^ , uniform perturbations in the rotated frame correspond to points along the line 
C 

A = c^fx/v , (24) 

in the (/x. A) plane. 

Whenever the evolution of an exponentially localized perturbation of type (^) is con- 
sidered, it is natural to introduce the quantity V = A//i, which can be interpreted as the 
velocity of the front . Eq. ( Plj) connects this velocity with that of the rotated frame, 

V = c^v . (25) 
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Therefore, on the basis of definition (p3D, {x',ct') can be interpreted as the reference frame 
in which the front associated with the perturbation propagates with an "infinite" velocity. 
The exphcit expression for rj is 

V = .JX' + (c/i)2 = X/(3 . (26) 

Such a relation can be turned into a self-consistent equation for the maximum Lyapunov 
exponent rj^ax by imposing the constraint that the pair (/z, A) lies on the line A = Amax(/^), 
namely 

'7max "^'^max ( ~^PVma,x | • (^'^) 



/5 "'^nc' 

Some ambiguities arise when velocities v > c are considered, since the line C intersecates 
Amax(/u) in two points as seen in Fig. ^. This phenomenon was already noticed in Ref. |]10 



while discussing the propagation of exponentially localized disturbances in the original ref- 
erence frame. Moreover, it has been shown that only the front corresponding to the smaller 
value of /i is stable, except for some cases where a nonlinear mechanism intervenes domi- 
nating the propagation process [11]. 

Atv = (? jV^ the two intersections degenerate into a single tangency point. This condition 
defines K, which can be interpreted as the slowest propagation velocity of initially localized 



disturbances 10 



The extension of Eq. (^) to the rest of the spectrum requires to connect n\ and with 
rZr,. In the next section, we will show how to perform such a step with the help of the entropy 
potential. Here, we limit ourselves to discuss the structure of the STLS for different values 
of the tilting angle In Fig. ^, we report the borders of the bands, which can be determined 
from the intersections of L with the border dT> of the domain of allowed perturbations (see 
Fig. 1 of LPT and Fig. ^). For -i? = (temporal case) a single band is present but, as soon 
as -i? > 0, a second negative band arises from the intersections with the branch diverging 
to — oo at /i = — /ic- For > 7r/4, the negative band disappears and a positive band arises 
from the intersections with the branch diverging to +oo with slope v = c. K single band 
spectrum is again recovered for ^9 > -i?* = atan(c/K)- 

Notice that in symplectic maps, the STLS is symmetric for any value of d (see LPT) so 
that positive and negative bands appear and disappear simultaneously. 

It is worthwhile to illustrate some of the above considerations in the simple case of the 
linear diffusion equation (^. The expression for Amax(Ai) can be obtained from Eq. (^, by 
setting k and uj equal to 0. Accordingly, Eq. ( pTj) reads as 

= 7 + ^ (^/^w) • (28) 

On the other hand, the model equation in the rotated frame can be obtained from the 
substitutions 

dt^fi{d,-vd^,) (29) 
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By introducing the usual Ansatz for the shape of the perturbation, 

^(x', t') ~ exp [mx' + {rj + iQ)t'] , (30) 

separating the real from the imaginary part, and eliminating Q, we obtain the integrated 
density of spatiotemporal exponents 

v)=/3(^l- 2D^(3v^ \/^(/3^)^-§ + ^ (31) 

and the corresponding STLS //(k, v). Dimensional analysis shows that n is an inverse length, 
as expected for a density of exponents. Notice that Eq. (^) is recovered, by setting k = 
n = Om Eq. (|31|). 

In this and more general continuous models, we should remark that the line C intersects 
Amax(Ai) twice for any arbitrarily small v. This is because the Laplacian operator sets no 
upper limit to the propagation velocity of disturbances. 

IV. FROM THE ENTROPY POTENTIAL TO DYNAMICAL INVARIANTS 

The present section is devoted to establish the consequences of the existence of the 
entropy potential on the Lyapunov spectra and other dynamical indicators such as the 
Kolmogorov- Sinai entropy and the Kaplan- Yorke dimension of the attractor. In order to 
keep the notations as simple as possible, we assume that time and space coordinates are 
scaled in such a way that c = 1. 

A. Spatiotemporal exponents 

The very existence of the entropy potential ^ implies that the Lyapunov spectrum in a 
frame tilted at an angle (recall that ^ is the angle from the A-axis) can be obtained by 
computing the derivative of $ along the direction u = (sint?, cos ??) in the (yU, A) plane. In 
fact, this is a straightforward generalization of the previous findings that ra^ and nx are the 
derivatives of $ along the fi and A direction, respectively. Accordingly, the STLS is linked 
to the TLS and SLS by the following general equation 

n^{v,r]) = u-V^ = f3[vnf, + nx] , (32) 

where V = {d^, d\) is the gradient in the (/x. A) plane, and the r.h.s of the above formula is 
evaluated for 

/i = vPt] (33) 
X = (3r] . (34) 

Such a relation can be directly verified for the diffusion equation from Eqs. (PD , (|3TD . Further, 
more significative tests have been performed by checking numerically the validity of Eq. (|32D 
in some lattice models involving, e.g., logistic and homogeneous chains, the spectra of which 
are reported in see Fig. |^ (notice that, since the time axis has not been renormalized in the 
tilted frame, the factor /3 need not be introduced). 
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B. Entropy 



Kolmogorov- Sinai entropy H^s is a measure of the information-production rate during 
a chaotic evolution. An estimate of Hks is given by the Pesin formula |jl2[ as the sum Hx of 
the positive Lyapunov exponents. While it is rigorously proven that Hks ^ H\, numerical 
simulations indicate that, in general, an equality holds. In spatially extended systems, H^s 
is believed to be proportional to the system size |T^. For this reason, it is convenient to 
introduce the entropy density hx which, in the thermodynamic limit, is computed as the 
integral of the positive part of the Lyapunov spectrum. 

Therefore, it is natural to extend the definition of entropy to tilted frames as an integral 
along the line £, 



h^= ri{n^,nx)dnjj , (35) 

where nmin is the integrated density in the point where the line C intersects &D, i.e. where 
V = Vmax, while n^ax is measured in the origin, i.e. where 77 = 0. In the limit v = 0, 
the above equation reduces to the definition of the density hx, which refers to the original 
reference frame. For v ^ 00, instead, we obtain the "spatial" entropy density /i^. 

Numerical simulations performed with different CML models indicate that hx < h^. 
This can be explained by the following argument. The patterns obtained asymptotically 
by iterating the model in the original reference frame are, in general, unstable if generated 
along the spatial direction [Q. In other words, the spatiotemporal attractor is a (strange) 
repellor of the spatial dynamics. Accordingly, part of the local instability accounted for by 
the sum of positive spatial Lyapunov exponents is turned into a contribution to the escape 
rate from the repellor and h^ must be larger than the entropy hx of the original pattern. 



In continuous models this inequality is brought to the extreme h^ is infinite. 

Numerical simulations of a lattice of logistic maps with < f < confirm the 

previous evidence that /i^ is independent of v. Moreover, it is found that /i^ is constant 
and equal to h^ for v large enough. Both results have a simple explanation in terms of the 
entropy potential. The very definition of hr^ suggests that it is more convenient to refer to 
the (n^, nx) plane. Indeed, as long as the outer intersection of C with dV in the plane (/i. A) 
occurs on the same branch, the endpoints of the integrals ( |5BD are the same (in fact, in Fig. 3 
of LPT it is shown that each connected component of dV is mapped onto a single point in 
the (n^, nx) plane). 

If we now notice that the integrand in Eq. ( pSj ) is the gradient of the potential 

<l = \nx + /in^ - $ , (36) 

it is mathematically obvious why /i^ is independent of v in suitable intervals. 

Another explanation of the above result can be found by referring directly to the plane 
(/i. A). Integrating by parts, Eq. (|35|) can be rewritten as 

hjj = nmin?7max - / ^.^(/i, A)rf?7 . (37) 



As ninin = along the upper border (see LPT), we can compute /i^(f ) for t> < 1 by integrating 
the gradient of $ along the line C from (0, 0) to the upper border itself. Since the upper 
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border is an equipotential line (see also Fig. |l]), /i^ is independent of v. The same argument 
can be repeated for v > 1/V* by suitably shifting and n^. 

The independency of /i^ of v has also a physical interpretation. The Kolmogorov-Sinai 
entropy density, in fact, is the amount of information needed to characterize a space-time pat- 
tern (apart from the information flow through the boundaries |13[) divided by its temporal 
duration and the spatial extension, i.e. divided by the area. Therefore, hxs is independent 
of they way the axes are oriented in the plane, i.e. of the velocity v. As a consequence, 
hfj = hxs for all f < 1. 

The above conclusion still holds when the STLS exhibits a positive band as well (which is 
always the case in continuous models), provided that the content of such a band is discarded. 
Accordingly, we can conclude that the new degrees of freedom, associated to the positive 
band, which appear in the rotated frame are just physically irrelevant directions which turn 
the original attractor into a repellor. If f > c^/K, the two bands merge together and it is not 
anymore possible to distinguish between unstable but irrelevant directions and the unstable 
manifold of the original attractor. Presumably, this means that the repellor is turned into 
a strange repellor with a singular measure along some (all) unstable directions. 



C. Dimension 



A second important indicator of the "complexity" of a spatiotemporal dynamics is the 



fractal dimension. An upper bound Dky to it is given by the Kaplan- Yorke formula |]T2 
The existence of a limit Lyapunov spectrum, implies that Dky is proportional to the system 
size |TB|, so that it is convenient to introduce the dimension density dxy- In the framework of 
the present paper, it is natural to extend the definition to generic velocities. The dimension 
density satisfies the integral equation, 

CY 

ri{v,nri)dnjj = . (38) 

As for the entropy density, Eq. (|38| ) can be more easily interpreted with reference to the 
{nfj_,n\) plane. In fact, the curve implicitely defined by the above constraint is the equipo- 
tential line C 

^{n^,nx) = . (39) 

The dimension density dxyiv) can, in turn, be determined from Eq. ( P^ at the intersection 
point between C and the image of C in the plane {n^,nx). 

At variance with the entropy density, dxyiv) changes with v (see Fig. P) even if we 
avoid considering the second positive band. In fact, while h^j is an information divided by 
a space-time area, dxriv) is a number of degrees of freedom divided by a length, measured 
orthogonally to the propagation axis. Thus, at least from a dimensional point of view, it is 
meaningless to compare dxriv) for different velocities. However, one can reduce temporal 
to spatial lengths by introducing the scaling factor c and, in turn, ask himself how the 
dimension changes with c. It is easily seen that the scaling dependence on c is expressed by 
the following relation. 
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dKY{v,Ci)^l+ i^ — j = dKY{v,C2)^l+ i^ — j , (40) 

The (completely arbitrary) choice of c reflects in different dependences of dxY on v. A 
natural procedure to fix c is by minimizing the dependence of dxv on the observation angle. 
This amounts to choosing the time units in such a way as to make the 2D pattern as isotropic 
as possible. In homogeneous CMLs, the procedure is so effective that a suitable choice of c 
allows removing alomost completely the velocity dependence as seen in Fig. |^a, where the 
results for the natural value c = 1 are compared with those for c = 3. 

More in general, however, it is not possible to achieve such a complete success. This is, 
for instance, the case of the logistic CML, where the dimension drop for c = 1 is too large 
to be compensated by any choice of c (see Fig. ^d, where the curve for c = 1 is compared 
with the best results obtained for c = +oo). 

A further indicator which is sometimes useful in characterizing the chaoticity of a given 
extended system is the dimension density du of the unstable manifold. This dimension is 
nothing but in the point where ?7 = 0, i.e. in the origin, and its expression simply reads 
as 

d^ = /3nxi0,0) . (41) 

being n^(0, 0) = 0. The choice c = +oo of the scaling factor removes exactly the dependence 
on the orientation of the reference frame. This choice is equivalent to measuring lengths in 
the untilted frame. 



V. COMOVING EXPONENTS 

Another class of indicators, introduced to describe convective instabilities in open-flow 
systems, consists of the so-called comoving or velocity-dependent Lyapunov exponents 0. 
They quantify the growth rate of a localized disturbance in a reference frame moving with 
constant velocity V. Given an initial perturbation u{x, 0) which is different from zero only 
within the spatial interval [—Lq/2, Lq/2], numerical analyses indicate 

M(x,t) ~exp(A(a;/t)t) , (42) 

for t sufficiently large. Eq. (^) defines the comoving Lyapunov exponent A as a function of 
V = x/t. The initial width Lq of the disturbance is not a relevant parameter, since a generic 
perturbation grows with the maximum rate.0 

The definition of A can be extended to a whole spectrum of comoving exponents by 
looking not just at the local amplitude of the perturbation but also at its shape |T^. Since 



the physical meaning of the rest of the spectrum is still questionable, in the following we 
limit ourselves to discuss the maximum. 



^In the particular case of a J-like initial profile, the definition of local Lyapunov exponent intro- 
duced in Ref. 11161 is recovered. 
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As a matter of fact, the limit t — ^ oo (required by a meaningful definition of an asymptotic 
rate) implies the infinite-size limit. Therefore, one must carefully keep under control the 
system size, when longer times are considered. This is perhaps the most severe limitation 
against an accurate direct measurement of A. 

It can be easily shown that A{V) is connected with the maximal temporal Lyapunov 
exponent Amax(/^) by a Legendre-type transformation fl^. Eq. (^2]) implies that the pertur- 
bation has a locally exponential profile with a rate 

dAjV) .... 

in the point x = Vt. On the other hand, we know that such a profile evolves as 

u{Vt, t) ~ exp[(A,„ax(/i) + I^V)t] . (44) 
By combining Eqs. (^21) and (0), we obtain 

A(y) = A^ax(/i)+/i^^^^^ , (45) 

which, together with Eq. (^) can be interpreted as a Legendre transform from the pair 
(A, V) to the pair (Amax, f^)- The inverse transform reveals the further constraint 

V = ^^Hldl^ . (46) 

Eq. ( P5| ) states that A(V") is the growth rate of an exponentially localized perturbation with 
a given fi value as determined from the condition Eq. (^31). However, the perturbation 
itself propagates with yet another velocity, V = Aniax(At)//U- As a matter of fact, V and V 
correspond to phase and group velocities for propagating waves in linear dispersive media. 
In particular, the "phase" velocity V can be larger than the "light" velocity (c = 1 in CML 
with nearest neighbour coupling), while V is bounded to be smaller. 

A simple geometrical interpretation of the above Legendre transformations can be given 
with reference to the (/z. A) plane. The comoving Lyapunov exponent A(V^) is the distance 
between the origin and the intersection of the A axis with the straight line of slope V, tangent 
to the upper temporal border. If the system is chaotic, such an intersection remains positive 
for V < V^:. Therefore, as already remarked, is the maximum velocity of disturbance 
propagation. Indeed, along the worldlines with V > V^,, the disturbance does not vanish 
exactly, but decreases exponentially in time. 

Whenever a Legendre transform comes into play, some attention must be payed to the 
concavity of the functions involved in the transformation. In the present context, this is 
the case of frozen random patterns where the border of the allowed region exhibits a change 
of concavity at /i = /ii (see Fig. 6 of LPT). This implies that for /i < /ii, the maximal 
temporal exponent is constant and equal to Aniax(Aii)- The corresponding "phase transition" 
reflects itself as a linear dependence of the comoving Lyapunov exponent on the velocity for 
\V\<V,, 

A{V) = X^^M - f^iV . (47) 
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This is evident in Fig. ^, where the whole set of A values is reported. 

It is clear that comoving, spatiotemporal and temporal exponents are related to one 
anothe. However, the link is not so straightforward as one might think. Indeed, the velocity 
V of the rotated frame where ?7max coincides (up to a normalization factor) with Amax is equal 
to /V and thus differs from both V and V . 

VI. CONCLUSIONS 

In the present paper we have shown that all instability properties of ID chaotic sys- 
tems can be derived from a suitable entropy potential expressed as a function of any pair 
of variables in the set {/x. A, n^, n^}- The most appropriate representation depends on the 
problem under investigation. For instance, the properties of Kolmogorov-Sinai entropy are 
more naturally described with reference to {n^^nx). This is analogous to standard thermo- 
dynamics, where several potentials (Gibbs, Helmholtz, etc.) are introduced to cope with 
different physical conditions. 

The very notion of entropy potential implies general relations among the classes of Lya- 
punov exponents introduced and discussed here and in LPT, namely spatial, temporal, 
spatiotemporal and comoving exponents. Another remarkable consequence of the existence 
of an entropy potential is the independency of the Kolmogorov-Sinai entropy density h^s 
(as determined from the the spatiotemporal spectrum) of the propagation direction in the 
space-time plane. Accordingly, hxs can be considered as a super-invariant dynamical in- 
dicator. This is not the case of the fractal dimension, the dependence of which provides 
information about the anisotropy of the pattern. 

We should, however, point out that our statements are not rigorously proved (except for 
some simple test models). However, since our numerical simulations suggest their general va- 
lidity, we strongly believe that systematic analytical investigations should eventually succeed 
in proving their validity. A final remark concerns the space dimensionality. The existence of 
the entropy potential stems from the analyticity of the complex dispersion relations which, 
in turn, is peculiar of ID systems. 
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APPENDIX A: 

The crucial point in justifying the existence of the entropy potential is the analytic 
structure of the eigenvalue equation stemming from the linearized dynamics. To support 
the generality of this statement we consider in this Appendix two more examples, namely 
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the linear stability analysis of homogeneous solutions both of the ID wave equation and of 
CML. 

The wave equation 

d^u = -rn^u + dlu, (Al) 

is the conservative analogous of Eq. (Q) (m is a real parameter) and can be treated in a 
similar way, obtaining 

The abovee expression justifies per se the existence of the entropy potential. Incidentally, 
notice that the Hamiltonian nature of Eq. ( [AID implies the degeneracy of the standard 
TLS in zero, since the uniform solution is an elliptic fixed point. The entropy potential is 
determined as the real or, equivalently, the imaginary part of the formal integral 



— cosh ^ f — 



m 



(A3) 



This can be verified in the limit of a "weak" instability m — > 0, when Eq. ([A3|) approximately 

reads as 



(A4) 



^(fi)^- /i^-m^logf-^ 
2 L Km, 

By also expanding to the lowest order in m the expressions of A and u determined by 
Eq. we obtain 

/ , N , / 1 \ 
uj{fi, k) ^ k [1 + -- 



2 ^2 + ^2 y 
It is straightforward to verify that 

dkRe-^ = -9^Im§ = u (A6) 
a^Re^ = dklm^ = X . 

For homogeneous solutions of CML models, we obtain 

= r [(1 -e) +£:cosh/i] , (A7) 

where r is the multiplier. Unfortunately, in this case it is not possible to write down an 
explicit expression for the integral for generic parameter values. We limit ourselves to 
discuss the problem in the limit of a small coupling, i.e. e ^ 0. Expansion of ( [A^ ) to the 
first order in e, yields 

^(/2) (logr — £:)/i + esinh/i , (A8) 

and 

A(/i, nx) ~ logr — e {1 — cos k cosh fi) (A9) 
n^{fi,nx) ^ £ sink sinh fi , 

which should be compared with the corresponding expressions obtained by expanding to first 
order in e Eqs. (16) and (20) of LPT. Moreover, one can verify that the relations analogous 
to Eqs. ( |A^ ) hold also in the present example. 
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FIGURES 



FIG. 1. Contour plot of the entropy potential ^> for a homogeneous chain (r = 2, e = 1/3). 

FIG. 2. Lattice implementation of the definition of spatiotemporal Lyapunov exponents for 
V = 1/2. 

FIG. 3. Plot of the boundary dT> and of the line A = u/i in the three velocity regimes for the 
logistic CML e = 1/3. The three lines refer to the different cases f < 1 (solid), 1 < v < 1/1^ 
(dot-dashed) and v = (dashed). 

FIG. 4. Boundaries of the STLS versus the tilting angle 'd for the homogeneous chain (r = 2, 
e = l/3). 

FIG. 5. Comparison between the STLS obtained by direct numerical computation and formula 
(|3^ ) for (a) homogeneous (r = 2) with w = 4/5 and (b) logistic CML with v = 3/5 (in both cases 

£ = 1/3). 

FIG. 6. Kaplan- Yorke dimension density cIky obtained from the STLS versus the tilting angle 
■!? for (a) homogeneous (r = 1.2) and (b) logistic CML models: in both cases e = 1/3. Circles refer 
to the scaling factor c = 1, while crosses correspond to c = 3, -|-oo in (a), (b), respectively. 

FIG. 7. Maximum comoving Lyapunov exponent A(y) for a frozen random pattern obtained 
as Legendre transform versus V (for comparison see also Fig. 6 in LPT). The vertical line indicates 
the position of the critical velocity Vi (see the text for definition). 



18 



TABLES 



Palli IiiLc-gral PmlIi IiiLc-gral <I> 

(0, 0) ^ (0, 3) 0.1011 (0, 3) ^ (2, 3) -0.3883 -0.2872 

(0, 0) ^ (2, 0) -0.5274 (2, 0) ^ (2, 3) 0.2406 -0.2868 

(0,0)^(0,0.8) 0.1069 (0,0.8)^(3,0.8) -0.8036 -0.6967 

(0,0)^(3,0) -1.4972 (3,0)^(3,0.8) 0.8 -0.6972 



TABLE I. Entropy potential ^> computed by integrating along two different paths in two 
different points of the (/x. A) plane. The difference is definitely smaller than the statistical error 
(« 10-3). 



19 



